Spectra and statistics of velocity and temperature fluctuations in turbulent 

convection. 
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Direct measurements of the velocity in turbulent convection of SFe near its gas-liquid critical 
point by light scattering on the critical density fluctuations were conducted. The temperature, 
velocity and cross frequency power spectra in a wide range of the Rayleigh and Prandtl numbers 
show scaling behaviour with indices, which are rather close to the Bolgiano-Obukhov scaling in 
the wavenumber domain. The statistics of the velocity fluctuations remains Gaussian up to the 
Reynolds numbers 10 . 



PACS numbers: 47.27.-i; 47.27.Te; 44.25.+f 

Turbulent convection exemplifies many startling as- 
pects of turbulent flows. The discovery of scaling laws in 
the heat transport and temperature statistics shed new 
light on the nature of the convective turbulence E10|. 
Bursts of thermal plumes from thermal boundary lay- 
ers and a coherent large scale circulation, which modi- 
fies the boundary layers via its shear, are found to co- 
exist in a closed convection cell [B-g|. These two salient 
features which reflect the competition between buoyancy 
and shear, were used in theoretical models to explain the 
observed scalings in the temperature field 0,H,||. The 
predictions which follow from these models for the ve- 
locity field, are still far from an experimental verification 

Therefore direct measurements of the velocity field would 
become a crucial step to test predictions of various the- 
oretical models, to understand the mechanism of energy 
vs. entropy transfer, to determine relative contribution 
of the thermal plumes and the large-scale circulation to 
the heat transport, and to study the mechanism of an 
onset and generation of the large-scale flow by intermit- 
tent thermal bursts. 

Lack of the velocity field measurements is explained by 
the fact that both conventional methods such as a hot- 
wire anemometry and a laser Doppler velocimetry(LDV), 
are not suitable for turbulent convection, as was com- 
monly agreed [pi. Indeed a hot-wire anemometry is use- 
less due to strong temperature fluctuations which accom- 
pany velocity fluctuations, and thus do not permit to 
correctly measure velocity. Fluctuations in a refractive 
index due to the large temperature fluctuations corrupt 
a signal and to drastically reduce a signal-to-noise ratio 
of LDV measurements because of wandering and defo- 
cusing of laser beams. Various attempts to measure the 
velocity field were rather limited due to severe experi- 
mental difficulties to directly measure the fluid velocity 
in convective flows in a wide dynamical range JqM . Novel 
direct velocity measurements by dual- and single-beam 



two-color cross-correlation spectroscopy technique || in 
rather limited range of the Rayleigh numbers, Ra, still do 
not provide information about the velocity spectra even 
in the frequency domain and direct information about 
statistics in a sufficiently wide dynamical range. 
Another problem, which is still controversial and requires 
further experimental test, is whether the scaling laws in 
the inertial range of the temperature and velocity power 
spectra in homogeneous turbulent convection are either 
those suggested by Kolmogorov, or those suggested by 
Bolgiano and Obukhov(BO) 0,0] ■ It was argued Q], that 
only the Kolmogorov scaling k^/ 3 for both the tempera- 
ture and velocity power spectra, where k is the wavenum- 
ber, is consistent. It means that the temperature is a 
passive scalar which just follows the velocity field. Other 
theories predict that depending on the parameters range 
either the Kolmogorov or BO scalings can be observed 
||. The latter gives a power law fc~ n / 5 for the veloc- 
ity fluctuations spectrum, and k~ 7 ' 5 for the temperature 
fluctuations spectrum. Here the temperature is not a pas- 
sive scalar but provides a potential energy on all scales. 
First measurements of the frequency spectra of the tem- 
perature fluctuations W and then many subsequent ex- 
periments H[| reveal the BO power law, although the 
inertial range was very limited, and the applicability of 
the Taylor hypothesis is not completely justified. The 
only measurements of the velocity spectra in a real space 
by a homodyne spectroscopy method in rather narrow 
range of Ra also confirm the BO predictions Q. How- 
ever, various numerical simulations and, particularly re- 



cent remarkable ones [10 support the Kolmogorov scal- 
ing and call for further experiments. 
Our experimental studies during the last several years 
convincingly demonstrated that the Rayleigh-Benard 
convection in a gas near the gas-liquid critical point (CP) 
is a very appropriate system to study high Ra turbulent 
convection [jlh,f2| ■ Singular behaviour of the thermody- 
namic and kinetic properties of the fluid near T c pro- 



vides the opportunity to reach both extremal values of 
the control parameter (Ra up to 5 • 10 14 has been reached 
in our experiment on SFq ) and to scan Prandtl num- 
ber, Pr, over an extremely wide range ( from one till 
several hundred). All these features make the system 
unique in this respect ITlUl3]. However, the most rel- 
evant aspect of the system to the results presented in 
this Letter, is our finding that the presence of the critical 
fluctuations provides us the possibility to conduct LDV 
measurements in a rather wide range of proximity to CP 
between 3 ■ 10 -4 and 10~ 2 in the reduced temperature 
t = (T — T c )/T c , where T is the mean cell temperature. 
It corresponds to variation in Pr from about 20 till 200. 
The upper limit in t is defined by a scattered light in- 
tensity. It is experimentally found that at r > 10~ 2 the 
photomultiplicr current is too low to recover the Dopplcr 
frequency shift. The lower limit is determined by mul- 
tiple scattering in a large convection cell as a result of 
which the light intensity at the detector is drastically 
reduced. The signal-to-noise ratio of the LDV measure- 
ments is increased with the closeness to CP until it drops 
rather abruptly at r < 3 ■ 10~ 4 |l]@. 
Critical density fluctuations moving with the flow, scatter 
light as effectively as seeding particles in a conventional 
LDV: the corresponding signal from a photomultiplier 
is of excellent quality Q. The light scattering occurs 
on the short-scale thermal density fluctuations advected 
by the large scale turbulent fluctuations. In the range 
of r used in the experiment, the scattering takes place 
mostly on the thermal density fluctuations. Indeed, due 
to smallness of the Mach number at maximum flow veloc- 
ities the density fluctuations due to turbulent flow were 
negligible. 

It is well known from conventional LDV measurements 
on seeding particles, that the presence of too many seed- 
ing particles in a scattering volume leads to significant 
corruption of the signal. So, the question is, how so large 
density of scatterers provides so clean signal. The light 
scattering on the critical fluctuations occurs differently 
from that on the seeding particles: the scattering am- 
plitude is random in both the amplitude absolute value 
and phase in space and time. The crucial point here 
is that the scattering occurs on sound phonons due to 
large fluid compressibility, rather than on individual fluc- 
tuations of the order of the correlation length. So the 
time variations of the pressure fluctuations are corre- 
lated even at great distances. In the range of r men- 
tioned above, the correlation length is less than 10~ 5 cm. 
On the other hand, the wavenumber of the light excited 
sound phonons is q = 2fcosin#/2, where 6 is the angle 
between the crosssed laser beams in the LDV system, 
and k is the light wavenumber. Then the corresponding 
sound wavelength is A = 2n/q w 3/im, and the frequency 
/ = qc s w 20 Mhz, where c s is the sound velocity. At 
this frequency and in this range of r the sound attenua- 
tion length will be by orders of magnitude larger than the 



correlation length JLJ , and turns out to be the length on 
which the coherent light scattering takes place. Since it is 
also less or of the order of the beam waist in the scattering 
volume, it means that at these conditions the scattering 
occurs just on a few coherent scattering regions like on 
a few particles. That is the reason for observation of a 
good quality signal. Of course, the LDV used the scat- 
tered signal rather close to the central Rayleigh line, at 
frequency differences much smaller than the correspond- 
ing Mandelshtam-Brillouin doublet frequency shift due 
to the phonon scattering Alo = f (several kHz compared 
with 20 MHz) E3. But nonlinear coupling between the 
amplitude and the phase of the scattering light produces 
sufficient signal at low frequency shifts to be observable. 
In this Letter we studied Pr dependence of spectra and 
statistical properties of the vertical velocity fluctuations 
at one location in the frequency domain using this novel 
LDV technique. In addition we also obtained the tem- 
perature fluctuations spectra in a wide range of Ra and 
Pr. 

The experiment we present here, was done with a high 
purity gas SF 6 (99.998%) in the vicinity of T c and at 
the critical density (p c — 730kg /m 3 ). This fluid was 
chosen due to the relatively low critical temperature 
(T c = 318.73 K) and pressure (P c = 37.7 bar) and well- 
known thermodynamic and kinetic properties far away 
and in the vicinity of CP. The set-up was described 
elsewhere JL1||12| ]. The cell is a box of a cross-section 
76x76 mrrriormed by 4 mm plexiglass walls, which are 
sandwiched between a Ni-plated mirror-polished copper 
bottom plate and a 19 mm thick sapphire top plate of 
L = 105 mm apart. 

The cell is placed inside the pressure vessel with two side 
thick plastic windows to withstand the pressure differ- 
ence up to 100 bar. So the cell had optical accesses from 
above through the sapphire window and from the sides. 
They were used for both shadowgraph flow visualization 
and for LDV. The pressure vessel was placed inside a 
water bath which was stabilized with rms temperature 
fluctuations at the level of 0.4 mK. The gas pressure 
was continuously measured with 1 mbar resolution by 
the absolute pressure gauge. Together with calibrated 
100ft platinum resistor thermometer they provide us the 
thermodynamic scale to define the critical parameters of 
the fluid and then to use the parametric equation of state 
developed recently for SFq |15[ . 

Local temperature measurements in a gas were made by 
three 125/im thermistors suspended on glass fibers in the 
interior of the cell (one at the center, and two about half 
way from the wall). Local vertical component velocity 
measurements at about L/4 from the bottom plate were 
conducted by using LDV on the critical density fluctua- 
tions 0@. 

A typical time series of the vertical velocity component 
of turbulent convection at Pr = 93 and the power spec- 
tra at Ra from 3 • 10 11 to 8 • 10 13 are shown in Fig.l. The 



spectra for other values of Pr look similar. The spectra 
at higher Ra are characterized by well-defined low fre- 
quency peaks related to the large scale circulation, by 
surprisingly large inertial range, and by relatively high 
level of white noise which significantly limits the dynam- 
ical range of the spectra. The noise level does not de- 
pend on the flow velocity and remains at the same level 
regardless of the heating power which drives the convec- 
tion flow, but does depend on r. So the signal-to-noise 
ratio increases with approaching the CP and with a ve- 
locity increase. As was shown, the source of the white 
noise is the phase noise due to the critical density fluc- 
tuations. The results of these studies will be published 
elsewhere |Lj]| . 

The random phase of the scattering signal causes an ad- 
ditional noise in the velocity measurements, and thus 
limits the obtaining the reliable data. Indeed, the mea- 
sured Doppler frequency shift can be presented as <j) — 

V ■ k + (j) t h , where V ■ k is the velocity component along 
the wavevector k, and 4>th is the frequency shift due to 
the thermal phase noise. Then the corelation function of 
the measured frequency shift, averaged over long time, 
is (#) = {(V ■ fc) 2 ) + (Mth) + 2{(V ■ k ■ <£ th ). This 
expression at adiabatically slow temporal variations of 

V ■ k reduces to (</></>) = ((V ■ fc) 2 ) + {(/>th4 l th), that leads 
to the same additive relation in the corresponding fre- 
quency spectrum. So, at V ■ k S> <ftth one reproduces 
the velocity from the measured Doppler frequency shift. 
Besides, the smallest energy containing eddies should be 
larger than both the correlation length and beam waist. 
If these conditions are satisfied, one can expect that the 
frequency shift of the scattered light is equal to the fluid 
velocity in the wavenumber units. 

From the fit in the inertial range (Fig.l) one finds well- 
defined power law behaviour. In fact the LDV measure- 
ments do not include the whole inertial range due to the 
noise limitation. So one cannot reach the dissipation re- 
gion. The scaling index obtained by the linear fit, similar 
to that shown in Fig.l, and averaged over all values of 
Ra and Pr, gives — 2.4 ± .2, which rather close to the BO 
scaling in the wavenumber domain |8[. 
From the velocity time series the probability density func- 
tion (PDF) of the velocity fluctuations normalized by 
the rms velocity, was also constructed. The PDF of 
the velocity fluctuations as well as the velocity differ- 
ences at a finite time lag were Gaussian (Fig. 2) for all 
Ra and Pr under studies. We investigated also the de- 
pendence of the rms vertical velocity fluctuations and 
mean vertical velocity on Ra and Pr. The latter scaling 
behaviour was reported elsewhere |0,0]. The former 
one is presented in Fig. 3 in Re(— V s L/v) vs Ra coordi- 
nates. Here V s is the rms of vertical velocity fluctuations, 
and v is the kinematic viscosity. The fit to the data for 
three values of Pr = 27,45, and 93 suggest the scaling 
Re = 0.005i?a 043±02 with no Pr dependence. The data 
for Pr = 190 deviate from this scaling behaviour and can 



be better fitted with the scaling exponent 0.34. The for- 
mer scaling exponent of Re vs Ra for lower Pr is rather 
close to the value obtained for water(Pr = 7) and for 
the theoretical prediction 3/7 of both the mixing zone Q| 
and the large scale flow theories Q. We would like to 
emphasize that Pr dependence of Re in the bulk is very 
different from that found for the large scale circulation 

The power spectra of the temperature fluctuations were 
measured both at the cell center and on the side( L/4 
from the wall) in a wide range of Ra and Pr. We used 
also the data which we took in the compressed gas SFq 
far away from CP at P = 20bar, T = 303X and the 
density p — 0.18<?/cm 3 , that corresponds to Pr = 0.9, 
and at P = 50bar, T = 323X and p = 1.07g/cm 3 , that 
corresponds to Pr = 1.5. The data far away from CP 
cover the range of Ra between 10 9 and 5 • 10 12 . Most of 
the power spectra did not show any reasonable range of 
a scaling law. At Pr above 50 and the highest Ra the 
power law behaviour became noticable. At Pr > 90 and 
the highest values of Ra the scaling region of about one 
decade in frequency was observed (Fig. 4). The scaling 
index was obtained by fitting the spectra by H 



P{.f) = 



Bf s 
Bf s (f/f c 



Mf/fo) 



if / < fc 
if / > fc 



where s, / c , and a are the fit parameters. The scaling 
index obtained by the fit and averaged over all spectra, 
is s = —1.45 ± 0.1, that is close to the BO power law for 
the temperature fluctuations spectra in the wavenumber 
presentation. 

Thus the measured scaling indices for both the velocity 
and temperature fluctuations in the frequency domain 
are rather close to the theoretically predicted BO scal- 
ing in the wavenumber presentation. As commonly ac- 
cepted, the Taylor hypothesis can be used to pass from 
the wavenumber to frequency domain. Since we mea- 
sured both the velocity and temperature spectra in the 
rcqion with non-zero mean velocity, the application of 
the Taylor hypothesis is better justified. From the veloc- 
ity and temperature time seria taken at the same spatial 
location, we constructed cross-correlator which in the 
Fourier frequency domain also showed scaling behaviour 
with the scaling index —1.85 ±0.1, that is close to the 
theoretically predicted value in the BO region Hj. 
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FIG. 1. Power spectra of velocity fluctuations for Ra from 
3 ■ 10 11 to 8 ■ 10 13 at Pr = 93. Solid line is a power law fit. 
The insert: Vertical velocity time series at Ra — 3 ■ 10 12 . 



FIG. 2. PDF of the velocity time series for Ra from 3 ■ 10 11 
to 8 • 10 13 at Pr = 93. Each PDF is normalized by the rms 
of each time series. The thick line is a Gaussian PDF. 



FIG. 3. Re for the rms vertical velocity fluctuations 
as a function of Ra for three values of Pr: circles-27, 
up-triangles-45, down-triangles-93. Solid line is a power law 
fit. 



FIG. 4. Power spectra of temperature fluctuations at the 
cell center at Pr = 93, Ra = 8.4 ■ 10 13 -top; Pr = 190, 



Ra 



1.5 • 10 -middle; Pr = 300, Ra = 3 • 10 -bottom. 



In each plot a line of the best fit by the model described in 
the text, is shown. 
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